AD-A257  887 


NASA  Contractor  Report  189717 
ICASE  Report  No.  92-52 


ICASE 


Spric 

electe J 

issej 

c  > 


0 


A  GENERAL  INTRODUCTION  TO  AEROACOUSTICS 
AND  ATMOSPHERIC  SOUND 


James  Lighthill 


Appio'Tod  tew  pnicac  fstetis*; 
0w)3ibuaoo  UpLlisuTsd 


Contract  No.  NAS  1-18605 
October  1992 

Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  Virginia  23681-0001 

Operated  by  the  Universities  Space  Research  Association 


IVIASA 

National  Aeronautics  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Virginia  23665-5225 


92-30806 

■lililii 


for 

I  VTIS  ~«M4i 
!  WfC  tAI^ 
^kMi'Uiuuj«a4 

A  GENERAL  INTRODUCTION  TO  AEROACOUSTICS - 

AND  ATMOSPHERIC  SOUND 


“-5^ 


James  Lighthill^ 

Department  of  Mathematics 
University  College  London 
Gower  Street,  London  WCIE  6BT 
UNITED  KINGDOM 


Olstrilwittan/ 
Availability  Co* 
tA»all  aad/or 
I  Special 


Dlst 


9^\\  I 


ABSTRACT 

This  paper  uses  a  single  unifying  principle  (based  upon  the  nonlinear  “momentum-flux” 
effects  produced  when  different  components  of  a  motion  transport  different  components  of 
its  momentum)  to  give  a  broad  scientific  background  to  several  aspects  of  the  interaction 
between  airflows  and  atmospheric  sound.  First,  it  treats  the  generation  of  sound  by  airflows 
of  many  different  types.  These  include,  for  example,  jet-like  flows  involving  converted  turbu¬ 
lent  motions  -  with  the  resulting  aeroacoustic  radiation  sensitively  dependent  on  the  Mach 
number  of  convection  -  and  they  include,  as  an  extreme  case,  the  supersonic  “boom”  (shock 
waves  generated  by  a  supersonically  converted  flow  pattern).  Next,  the  paper  analyses  sound 
propagation  through  nonuniformly  moving  airflows,  and  quantifies  the  exchange  of  energy 
between  flow  and  sound;  while,  finally,  it  turns  to  problems  of  how  sound  waves  “on  their 
own”  may  generate  the  airflows  known  as  acoustic  streaming. 


'This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No,  NASl-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681-0001.  Research 
was  also  supported  by  the  Leverhulme  Trust. 


1.  Broad  Overview 


This  general  introductory  paper  is  devoted  to  Interactions  of  Sound  with  Air,  including 
transmission  through  the  atmosphere  and 


both  generation  of  sound  by 
and  propagation  of  sound  in 


I  airflows 


(e.g.,  manmade  flows  -  around  aircraft  or  air  machinery  -  or  natural  winds)  as  affected  by 
the  air’s  boundaries  and  atmospheric  composition;  with  (conversely)  generation  of  airflows 
by  sound  (acoustic  streaming). 

From  linear  acoustics  I  utilize  the  properties  of  the  wave  equation,  including 


(i)  the  short- wavelength  ray-acoustics  approximation  (Lighthill,  1978b,  hereinafter  de¬ 
noted  WF,  pp.  67-)  and 


(ii)  the  theory  of  multipole  sources  (WF,  pp.  31-)  -  with  the  long-wavelength  “compact- 
source”  approximation  (source  region  of  size  ^  with  uJ^|c  small,  where  u;  =  radian 
frequency,  radiates  like  a  concentrated  source); 

while  from  nonlinear  acoustics  I  use  (WF,  pp.  150-)  the  physics  of  waveform  shearing  and 
shock  formation. 

Techniques  special  to  aeroacoustics  and  atmospheric  sound  are  centered  on  the  momen¬ 
tum  equation  for  air.  Its  difference  from  a  wave-equation  approximation  include 

A.  Linear  effects,  of  gravity  acting  on  air  stratified  as  meteorologists  observe;  effects  which 
allow  independent  propagation  (WF,  pp.  292-)  of  “internal”  gravity  waves  and  of 
sound,  except  at  wavelengths  of  many  kilometers  when  the  atmosphere  becomes  a 
waveguide  (WF,  pp.  425-)  for  global  propagation  of  interactive  acoustic-gravity  waves; 

and  (still  more  importantly)  include 

B.  Nonlinear  effects,  of  the  momentum  flux  pUiUj;  i.e.  the  flux  -  rate  of  transport  across 
unit  area  -  of  any  pu,  momentum  component  by  any  Uj  velocity  component.  This 
term,  neglected  in  linear  acoustics,  acts  like  a  stress  (i.e.  force  per  unit  area  -  since 
rate  of  change  of  momentum  is  force).  In  particular. 


(i)  an  airflow’s  momentum  flux  puiUj  generates  sound  like  a  distribution  of  (time- 
varying)  imposed  stresses;  thus  not  only  do  forces  between  the  airflow  and  its 
boundary  radiate  sound  as  distributed  dipoles,  (— >  F  ~  — t-  F)  but  also  such 
stresses  (acting  on  fluid  elements  with  equal  and  opposite  dipole-like  forces)  ra¬ 
diate  (WF,  pp.  63-)  as  distributed  quadrupoles;  (— ><—  ~  — |-  4-  — )  (Lighthill, 
1952  and  Lighthill,  1962) 


(ii)  the  mean  momentum  flux  <  pUiUj  >  in  any  sound  waves  propagating  through 
a  sheared  flow  (with  shear  dVildxj)  is  a  stress  on  that  flow  (Lighthill,  1972  and 
WF,  pp.  329-),  and  the  consequent  energy  exchange  (from  sound  to  flow  when 
positive,  vice  versa  when  negative)  is 

<  puiUj  >  dVi/dxj;  (1) 

(iii)  even  without  any  pre-existing  flow',  energy-flux  attenuation  in  a  sound  wave  al¬ 
lows  streaming  to  he  generated  (WF,  pp.  337-)  by  unbalanced  stresses  due  to 
a  corresponding  attenuation  in  acoustic  momentum  flux  -  essentially,  then,  as 
acoustic  energy  flux  is  dissipated  into  heat,  any  associated  acoustic  momentum 
flux  is  transformed  into  a  mean  motion  (Lighthill,  1978a  and  WF). 

And  another  (less  crucial)  momentum-equation/ wave-equation  difference  is 

C.  Nonlinear  deviation  of  pressure  excess  p  —  po  from  a  constant  multiple,  Cq(p  —  po),  of 

density  excess. 

(i)  For  sound  generation  by  airflows,  this  adds  an  isentropic  term  to  the  quadrupole 
strength  per  unit  volume  (WF,  pp.  60-) 

Tij  =  puiUj  -1-  [(p  -  Po)  “  4ip  -  /’o)]  Sij,  (2) 

the  last  term  being  considered  important  mainly  for  flows  at  above-ambient  tem¬ 
peratures  (Lighthill,  1952  and 

Lighthill,  1963); 

(ii)  for  propagation  of  sound  with  energy  density  E  through  flows  of  air  with  adiabatic 
index  7,  the  mean  deviation  is  about  (7  —  l)E,  and  the  total  radiation  stress 
(Bretherton  and  Garrett,  1968) 

<  puiUj  >  +^(7  -  i)ESij  (3) 

adds  an  isotropic  pressure  excess  to  the  mean  momentum  flux  (although  the 
energy  exchange  (1)  is  unchanged  in  typical  cases  with  dVifdxi  essentially  zero). 

(And  we  note  that  the  very  special  case  of  sound  waves  interacting  on  themselves  [in  other 
words,  nonlinear  acoustics]  may  be  interpreted  (WF,  p.  148)  as  a  combined  operation  of  the 
“self-convection  effect”  B  and  the  (smaller)  “sound-speed  deviation”  C.) 
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2.  Compact  Source  Regions 

2.1.  Sound  generation  by  low-Mach-number  airflows 

The  main  nondimensional  parameters  governing  airflows  of  characteristic  speed  U  and 
length-scale  L  are  Mach  number  M  =  U/c  (where,  in  aeroacoustics,  c  is  taken  as  the  sound 
speed  in  the  atmosphere  into  which  sound  radiates)  and  Reynolds  number  R  =  U Ljv  where 
V  =  kinematic  viscosity.  Low-Mach-number  airflows  are  compact  sources  of  sound,  with 
frequencies 

regular  flow 
oscillations; 
extremely  irregular 
turbulence. 

Since,  in  either  case,  a  typical  frequency  uj  scales  as  Uf  L  (Strouhal  scaling),  the  compactness 
condition  ujL/c  small  is  satisfied  if  M  =  U/c  is  small  (Lighthill,  1962). 

A  solid  body  which,  because  of  flow  instability,  is  subjected  to  a  fluctuating  aerodynamic 
force  F  scaling  as  pU^L^  (at  frequencies  scaling  as  U/L),  radiates  as  an  acoustic  dipole  of 
strength  F,  with  mean  radiated  power  <  F^  >  /Vl-Kpi?. 

This  acoustic  power  scales  as  pU^l?  f<?  (a  sixth-power  dependence  on  flow  speed).  There¬ 
fore  acoustic  efficiency,  defined  as  the  ratio  of  acoustic  power  to  a  rate  of  delivery  (scaling 
as  pU^l?)  of  energy  to  the  flow,  scales  as  (U/c)^  =  M^. 

(Exceptions  to  compactness  include  bodies  of  high  aspect-ratio;  thus,  a  long  wire  in 
a  wind  [where  the  scale  L  determining  frequency  is  its  diameter]  radiates  as  a  lengthwise 
distribution  of  dipoles.) 

Away  from  any  solid  body  a  compact  flow  (oscillating  or  turbulent,  with  frequencies 
scaling  as  U/L)  leads  to  quadrupole  radiation  (see  B(i)  above)  with  total  quadrupole  strength 
scaling  as  pU^L^.  Acoustic  power  then  scales  as  pU^L^/c^:  an  eighth-power  dependence 
(Lighthill,  1952  and  Lighthill,  1962)  on  flow  speed.  In  this  case  acoustic  efficiency  (see 
above)  scales  as  (U/c)^  =  A/®. 

Such  quadrupole  radiation,  though  often  important,  may  become  negligible  near  a  solid 
body  when  dipole  radiation  due  to  fluctuating  body  force  (with  its  sixth-power  dependence) 
is  also  present  (Curie,  1955  and  Ffowcs  Williams  and  Hawkings,  1969). 

Near  not  necessarily  compact  bodies  a  more  refined  calculation  -  using  Green’s  functions 
not  for  free  space  but  for  internally  bounded  space  -  leads  in  general  to  the  same  conclu¬ 
sion;  that  quadrupole  radiation  with  its  eighth-power  dependence  is  negligible  alongside  the 
sixth-power  dependence  of  dipole  radiation  due  to  fluctuating  body  forces;  but  important 
exceptions  to  this  rule  include  sharp-edged  bodies,  where  features  of  the  relevant  Green’s 
function  imply  a  fifth-power  dependence  on  flow  speed  of  acoustic  radiation  from  turbulence 


narrow-banded 
at  moderate  R 
broad-banded 
at  high  R 


when  flow 
instabilities 
lead  to 


.3 


(Ffowcs  Williams  and  Hall,  1970;  Crighton  and  Leppington,  1970;  Crighton  and  Leppington, 
1971;  and  Crighton,  1981). 

2.2.  Sound  generation  by  turbulence  at  not  so  low  Mach  number 

The  chaotic  character  of  turbulent  flow  fields  implies  that  velocity  fluctuations  at  points 
P  and  Q,  although  they  are  well  correlated  when  P  and  Q  are  very  close,  become  almost 
uncorrelated  when  P  and  Q  are  not  close  to  one  another. 

Reminder:  statisticians  define  correlation  coefficient  C  for  the  velocities  up  and  Uq  as 
C  =<  vp  ■  vq  >  f  <  Vp  Vq  in  terms  of  the  deviations,  vp  =  up—  <  up  > 

and  Vq  =  Uq—  <  Uq  >,  from  their  means.  When  two  uncorrelated  quantities  are  combined, 
their  mean  square  deviations  are  added  up: 

<  Vp  -f  Vq  >^=<  Vp  >  +2  <  Vp  •  Vq  >  -f  <  Vq  > 

=<  Vp  >  +  <  Vq  >  if  (7  =  0. 

Theories  of  turbulence  define  a  correlation  length  with 

,  f  well  correlated  1  ~  ,  f  1  1 

up  a.nd  u«  I  I  -  -C  close  j  q  /  “  ' 

when  PQ  is  substantially 

Roughly  speaking,  different  regions  of  size  (  (“eddies”)  generate  sound  independently,  and 
the  mean  square  radiated  noise  is  the  sum  of  the  mean  square  outputs  from  all  the  regions 
(Lighthill,  1954).  Typical  frequencies  in  the  turbulence  are  of  order  u)  =  v/i,  where  v 
is  a  typical  root  mean  square  velocity  deviation  <  v^  so  that  for  each  region  the 

compactness  condition  ui/c  small  is  satisfied  if  v/c  is  small.  Compactness,  then,  requires 
only  that  a  r.m.s.  velocity  deviation  v  (rather  than  a  characteristic  mean  velocity  U)  be 
small  compared  with  c  -  which  is  less  of  a  restriction  on  M  =  U fc  and  can  be  satisfied  at 
“not  so  low”  Mach  number  (Lighthill,  1963). 

3.  Doppler  Effect 

How  is  the  radiation  from  such  “eddies”  modified  by  the  fact  that  they  are  being  converted 
at  “not  so  low”  Mach  number?  The  expression  Doppler  effect,  covering  all  aspects  of  how 
the  movement  of  sources  of  sound  alters  their  radiation  patterns,  comprises  (i)  frequency 
changes  (WF),  (ii)  volume  changes  (Lighthill,  1952  and  Lighthill,  1962),  (iii)  compactness 
changes  (Lighthill,  1963  and  Ffowcs  Williams,  1963). 
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3.1.  Frequency  changes 

When  a  source  of  sound  at  frequency  a;  approaches  an  observer  at  velocity  tu,  then  in  a 
single  period  T  =  2ic I ui  sound  emitted  at  the  beginning  travels  a  distance  cT  while  at  the 
end  of  the  period  sound  is  being  emitted  from  a  source  that  is  closer  by  a  distance  wT.  The 
wavelength  A  (distance  between  crests)  is  reduced  to 

\  —  cT  —  wT  =  27r(c  —  tn)/u;  (4) 

(Figure  1)  and  the  frequency  heard  by  the  observer  (27r  divided  by  the  time  A/c  between 
arrival  of  crests)  is  increased  to  the  Doppler-shifted  value  (WF) 

(jJ 

uJr  =  - z — r  :  the  relative  frequency  (5) 

1  —  [wfc) 

that  results  from  relative  motion  between  source  and  observer.  For  an  observer  located  on  a 
line  making  an  angle  9  with  a  source’s  direction  of  motion  at  speed  F,  the  source’s  velocity  of 
approach  towards  the  observer  is  lo  =  F  cos^  (Figure  1)  and  the  relative  frequency  becomes 

w  f  augmented  1  ,  .  .  f  acute  1  ,  ... 

=  l-(V7c)cos«  ^  1  diminished  )  ^  |  obtuse  / 

Such  Doppler  shifts  in  frequency  are  familiar  everyday  experiences. 

3.2.  Volume  changes 

When  an  observer  is  approached  at  velocity  u;  by  a  source  whose  dimension  (in  the 
direction  of  the  observer)  is  sounds  arriving  simultaneously  (Figure  2)  from  the  source’s 

I  imar  |  have  been  emitted  |  |  by  a  time  r  (say). 

In  the  time  t  for  sound  from  the  far  side  to  reach  the  observer,  after  travelling  a  distance 
ct,  the  relative  distance  of  the  near  side  in  the  direction  of  the  observer  was  increased  from 
^  to  ^  -f-  cur  before  it  emitted  sound  which  then  travelled  a  distance  c{t  —  r).  Both  sounds 
arrive  simultaneously  if 

ct  =  ^  -\-  WT  -f  c{t  —  r),  giving  r  =  - 


and  £  +  WT  =  - ; — =  iwr/u). 

1  —  [wfc) 


The  source’s  effective  volume  during  emission  is  increased,  then,  by  the  Doppler  factor  cur/cu 
(since  dimension  in  the  direction  of  the  observer  is  so  increased  whilst  other  dimensions  are 
unaltered:  Lighthill,  1952  and  Lighthill,  1962). 


If  turbulent  “eddies”  are  effectively  being  convected,  relative  to  the  air  into  which  they 
are  radiating,  at  velocity  V,  then  Equation  (6)  gives,  for  radiation  at  angle  6,  the  Doppler 
factor  uJrloj  which  modifies  both  the  frequencies  at  which  they  radiate  and  the  effective 
volume  occupied  by  a  radiating  eddy. 

But  Equation  (2)  specifies  the  quadrupole  strength  T,j  per  unit  volume  for  such  an  eddy. 
Without  convection  the  pattern  of  acoustic  intensity  around  a  compact  eddy  of  volume  P 
and  quadrupole  strength  (^Tij  would  be 

<  {i^fijXiXjr~‘^)'^  >  /16ff^r^/?oc^;  (8) 

and,  since  different  eddies  of  volume  P  radiate  independently,  we  can  simply  add  up  mean 
squares  in  the  corresponding  expressions  for  their  far-held  intensities.  This  gives 

<  {fijXiXjr~'^Y  >  /IhTT^r^/JoC®  (9) 

as  the  intensity  pattern  radiated  by  unit  volume  of  turbulence.  The  Doppler  effect  modihes 
this,  when  the  compactness  condition  is  satished,  by  hve  factors  Ur/u)  (one  for  the  change 
in  source  volume  P  and  four  for  the  frequency  change  as  it  affects  the  mean  square  of  a 
multiple  of  the  second  time-derivative  of  Tij)  and  this  intensity  modihcation  by  a  factor 

[1 -(I//c)cos^]-5  (10) 

brings  about  an  important  preference  for  forward  emission  (Ffowcs  Williams,  1963). 

3.3.  Compactness  changes 

As  {V/c)  increases,  however,  the  Doppler  effect  tends  to  degrade  the  compactness  of 
aeroacoustic  sources  in  relation  to  forward  emission.  Not  only  does  w^/c  increase  in  propor¬ 
tion  to  Mach  number,  but  an  even  greater  value  is  taken  by  the  ratio  which  must  be 

small  if  convected  sources  are  to  be  compact.  A  restriction  on  the  extent  (10)  of  intensity 
enhancement  for  forward  emission  as  Vjc  increases  is  placed  by  these  tendencies  (Lighthill, 
1963;  Ffowcs  Williams,  1963;  and  Dowling  et  al.,  1978). 

They  can  develop,  indeed,  to  a  point  where  the  compact-source  approximation  may 
appropriately  be  replaced  by  its  opposite  extreme;  the  ray-acoustics  approximation  (WF). 
Thus,  for  supersonic  source  convection  {Yjo  1),  the  relative  frequency  (6)  becomes  infinite 
in 

the  Mach  direction  9  =  cos~’(c/K),  (II) 

and  radiation  from  the  source  proceeds  (Figure  3)  along  rays  emitted  at  this  angle  (Ffowcs 
Williams  and  Maidanik,  1965). 
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Explanatory  note;  the  source’s  velocity  of  approach  w  towards  an  observer  positioned 
at  an  angle  (1 1)  to  its  direction  of  motion  is  the  sound  speed  c;  thus,  not  only  is  the  generated 
wavelength  (4)  reduced  indefinitely  (the  ray-acoustics  limit)  but,  essentially,  different  parts 
of  a  signal  are  observed  simultaneously:  the  condition  of  stationary  phase  satisfied  on  rays 
(WF).  Sounds  emitted  (WF,  p.  196)  by  a  source  approaching  an  observer  at  a  speed  u> 
exceeding  c  are  heard  by  him  in  reverse  order  (“pap  pep  pip  pop  pup”  becomes  “pup  pop 
pip  pep  pap”!)  but  when  u;  =  c  all  the  sound  (vowels  and  consonants!)  are  heard  together 
as  one  single  “boom.” 

Further  note:  the  influences  placing  a  limit  on  the  signal  propagated  along  rays  may 
include  the  duration  S  of  well-correlated  emission  from  turbulent  “eddies”;  and,  also,  may 
include  nonlinear  effects  (see  §4.2.  Supersonic  booms). 


3.4.  Uniformly  valid  Doppler-efFect  approximations 


.Just  as  a  correlation  length  (.  for  turbulence  was  specified  in  §2.2,  so  a  correlation  duration 
S  can  be  characterized  by  the  requirement  that  moving  eddies  have  |  1 


velocities  at  times  differing  by  substantially 


Combined  use  of  correlation  length  ( 


and  duration  6  affords  an  approximation  to  the  radiation  pattern  from  convected  “eddies” 
that  has  some  value  at  all  Mach  numbers,  spanning  the  areas  of  applicability  of  the  compact- 
source  and  ray-acoustics  approximations. 

Figure  4  uses  space-time  diagrams  where  the  space-coordinate  (abscissa)  is  distance  in  the 
direction  of  the  observer.  Diagram  (a)  for  unconvected  “eddies”  approximates  the  region  of 
good  correlation  as  an  ellipse  with  axes  I  (in  the  space  direction)  and  ^  (in  the  time  direction). 
Diagram  (b)  shows  such  a  region  for  convected  “eddies”  whose  velocity  of  approach  towards 
the  observer  is  le;  thus,  it  is  Diagram  (a)  sheared  by  distance  w  per  unit  time. 

Signals  from  far  points  F  and  near  points  A^,  in  either  case,  reach  the  observer  simulta¬ 
neously  -  as  do  signals  from  other  points  on  the  line  F N  -  this  line  slopes  by  distance  c 
(the  sound  speed)  per  unit  time. 

Compact-source  case  (i)  with  w/c  small:  the  space  component  of  FN  in  Diagram  (b)  is 
^[1  —  (re/c)]-’,  just  as  in  Equation  (7)  for  normal  Doppler  effect  (neglecting  finite  6). 

Ray-acoustics  case  (ii)  with  w/c  =  I:  the  space  component  of  FN  is  cS. 

Intermediate  case  (iii)  with  lu/c  “moderately”  <  1:  the  space  component  of  FN  is  i 
multiplied  by  an  enhancement  factor 


(12) 


which  represents  the  effective  augmentation  of  source  volume  due  to  convection  (Ffowcs 


Williams,  1963). 

This  enhancement  factor  (T2)  is  applied  not  only  to  the  volume  term  in  the  cpiadrupole 
field  (9)  but  also  twice  to  each  of  the  pair  of  twice-differentiated  terms  inside  the  mean  square; 
essentially,  because  time-differentiations  in  quadrupole  fields  arise  (W’F)  from  differences 
in  the  time  of  emission  by  different  parts  of  the  quadrupole  source  region  (and  the  time 
component  of  FN  in  Diagram  (b)  is  simply  the  space  component  divided  by  c).  As  before, 
then,  five  separate  factors  (12)  enhance  the  intensity  field:  and,  with  w  replaced  by  I'cosf?, 
expression  (10)  for  the  overall  intensity  modification  factor  is  replaced  by 

{[l-(V7c)cost?p  +  (f/c^)'^}-^/^  (13) 

riiis  modification  factor  (13)  affords  us  an  improved  description  of  the  influence  of  Doppler 
effect  not  only  on  the  preference  for  forw’ard  emission  but  also  on  the  overall  acoustic  power 
output  from  converted  turbulence  (Lighthill,  1963  and  Ffowcs  Williams,  1963).  For  example. 
Diagram  (c)  gives  (plain  line)  a  log-log  plot  of  the  average  (spherical  mean)  of  (13)  as  a 
function  of  V/c  on  the  reasonable  assumption  that  (  =  0.6V^(5.  As  V/r  increases  this  average 
modification  factor  rises  a  little  at  first,  but  falls  drastically  like  .'3(V'/c)“''’  for  K/c  significantly 
>  1. 

Now  low-Mach-number  turbulence  away  from  solid  boundaries  (§2.1)  should  radiate 
sound  with  an  acoustic  efficiency  scaling  as  (U/c)^  where  if  is  a  characteristic  velocity  in  the 
flow.  With  V  taken  as  that  characteristic  velocity  (although  in  a  jet  a  typical  velocity  V'  of 
ed(iy  convection  would  be  between  0..5  and  0.6  times  the  jet  exit  speed),  the  modification  of 
(say)  an  acoustic  efficiency  of  10“^(V7c)^  for  low  Mach  number  by  the  average  modification 
factor  would  cause  acoustic  efficiency  to  follow  the  broken-line  curve  in  Diagram  (c),  tending 
asymptotically  to  a  constant  value,  O.OO.'j,  (aeroacoustic  saturation)  at  high  Mach  number. 
Such  a  tendency  is  often  observed  for  sound  radiation  from  “properly  expanded”  supersonic 
jets  (see  below). 

4.  Introduction  to  Aircraft  Noise 
4.1.  Aero-engine  and  airframe  noise 

How  are  aeroacoustic  principles  applied  to  practical  problems  -  such  as  those  of  studying 
aircraft  noi.se  with  a  view  to  its  reduction  (Crighton,  197.5;  Goldstein,  1976;  and  Goldstein, 
1984)? 

In  any  analysis  of  the  generation  of  sound  by  airflows,  we  may  need  first  of  all  to  ask 
whether  the  geometry  of  the  problem  has  features  that  tend  to  promote  resonance.  For 
example,  a  long  wire  in  a  wind  (§2.1)  generates  most  sound  when  vortex-shedding  frequencies 
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{  are  fairly  close  to  1  wire's  lowest  natural  frequency  of  vibration;  giving  good 
and  so  can  lock  on  to  J  '  <500 

correlation  of  sideforces,  and  so  also  of  dipole  strengths,  all  along  the  wire. 

Again,  a  jet  emerging  from  a  thin  slit  may  interact  with  a  downstream  edge  (parallel 
to  the  slit)  in  a  resonant  way  (Curie,  1953  and  Powell,  19.53b);  with  very  small  directional 
disturbances  at  the  jet  orifice  being  amplified  by  flow  instability  as  they  move  downstream 
to  the  edge,  where  they  produce  angle-of-attack  variations.  Dipole  fields  a.s.sociated  with  the 
resulting  sideforces  can  at  particular  frequencies  renew  the  directional  disturbances  at  the 
orifice  with  the  right  phase  to  produce  a  resonant  oscillation.  Some  musical  wind  instruments 
utilize  such  jet-edge  resonance’s,  reinforced  by  coincidence  with  standing-wave  resonances  in 
an  adjacent  pipe. 

But  in  the  absence  of  such  resonances  (leading  to  enhanced  acoustic  generation  at  fairly 
well  defined  frequencies)  airflows  tend  to  generate  acoustic  “noise''  whose  reaction  on  the 
flow  instability  phenomena  themselves  is  negligible. 

Resonances  analogous  to  the  above  which  need  to  be  avoided  in  aircraft  design  include, 
for  example, 

(a)  panel  flutter,  generated  at  a  characteristic  frequency  as  an  unstable  vibration  of  a 
structural  panel  in  the  presence  of  an  adjacent  airflow  (Dowell,  1975); 

(b)  screeching  of  supersonic  jets  from  nozzles  which,  instead  of  being  “properly  expanded” 
so  that  an  essentially  parallel  jet  emerges,  produce  a  jet  in  an  initially  non-parallel  form 
followed  by  shock  waves  in  the  well  known  recurrent  “diamond”  shock-cell  pattern;  the 
first  of  these,  replacing  the  edge  in  the  above  description,  can  through  a  similar  feedback 
of  disturbances  to  the  jet  orifice  generate  a  powerful  resonant  oscillation  (Powell,  19.53a; 
Powell,  19.53b;  and  Howe  and  Ffowcs  Williams,  1978). 

And  undesirable  resonances  may  also  be  associated  with  aeroengine  combustion  processes 
(C’andel  and  Poinsot,  1988).  But  we  turn  now  to  the  aircraft  noise  of  a  broad-banded 
character  that  remains  even  when  resonances  have  been  avoided. 

Then  aero-engine  jet  noise  proper  (Lighthill,  1963)  (that  is,  the  part  unrelated  to  any 
interaction  of  jet  turbulence  with  .solid  boundaries)  tends  to  follow  a  broad  trend  similar  to 
that  in  Figure  4;  where,  however,  becau.se  the  eddy  convection  velocity  V  is  between  0.5 
and  0.6  times  the  jet  exit  speed  (/,  the  acoustic  efficiency  makes  a  transition  between  a 
value  of  around  10~‘’M'’  in  order-of-magnitude  terms  for  subsonic  values  of  M  =  U/c  and 
an  asymptotically  constant  value  of  10“^  or  a  little  le.ss  for  M  exceeding  about  2. 

The  above  tendency  for  M  <  1  implies  that  noi.se  emission  from  jet  engines  may  be 
greatly  diminished  if  a  given  engine  power  can  be  achieved  with  a  substantially  lower  jet 
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exit  speed,  requiring  of  course  a  correspondingly  larger  jet  diameter,  L.  Furthermore,  with 
acoustic  power  output  scaling  as  p[f^L^/c^  (§2-1)  and  jet  thrust  as  noise  emission  for 

given  thrust  can  be  greatly  reduced  if  i'  can  be  decreased  and  L  increased  by  comparable 
factors. 

dVends  (along  these  lines)  in  aero-engine  design  towards  large  turbofan  engines  with 
higher  and  higher  bypass  ratios,  generating  very  wide  jets  at  relatively  modest  mean  Mach 
numbers,  have  massively  contributed  to  jet  noise  suppression  (whilst  also  winning  advantages 
of  reduced  fuel  consumption).  On  the  other  hand,  such  successes  in  suppressing  jet  noi.se 
proper  (originally,  the  main  component  of  noise  from  jet  aircraft)  led  to  needs  for  a  dedicated 
focusing  of  attention  upon  parallel  reductions  of  other  aircraft-noise  sources  (Crighton.  1972): 

(a)  those  associated  with  the  interaction  of  jet  turbulence  with  solid  boundaries  -  where 
sharp-edged  boundaries  (§2.1)  pose  a  particular  threat; 

(b)  fan  noise  emerging  from  the  front  of  the  engine  and  turbine  noise  from  the  rear; 

(c)  airframe  noise  including  acoustic  radiation  from  boundary-layer  turbulence  and  from 
interaction  of  that  turbulence  with  aerodynamic  surfaces  for  control  purposes  or  lift 
enhancement. 

.Some  key  areas  of  modern  research  on  aero-engine  and  airframe  noise  are: 

for  jet  noise,  techniques  for  relating  acoustic  output  to  vorticity  distributions  (Powell, 
1964  and  Mohring,  1978),  and  to  any  coherent  structures  (Ribner,  1964  and  Ffowcs  Williams 
and  Kempton,  1978),  in  jet  turbulence;  and  for  taking  into  account  (cf  §5.4)  propagation 
through  the  sheared  flow  in  a  wide  jet  (Phillips,  1960  and  Mani,  1976); 

for  noi.se  from  fans  and  propellers,  mathematically  sophisticated  ways  of  reliably  esti¬ 
mating  the  extr  I  of  cancellation  of  dipole  radiation  from  different  parts  of  a  rotating-blade 
system  (alongside  a  good  independent  estimate  of  quadrupole  radiation:  Parry  and  Crighton, 
1991): 

for  airframe  noi.se,  a  recognition  (Powell,  1960  and  Crighton,  1984)  that  massive  cancel¬ 
lations  act  to  minimize  noise  radiation  from  bound  ary- layer  turbulence  on  a  flat  surface  of 
uniform  compliance  -  and,  therefore,  that  avo'dance  of  sharp  nonuniformities  in  airframe 
skin  compliance  may  promote  noise  reduction. 

4.2.  Supersonic  booms 

In  addition  to  aero-engine  and  airframe  noise,  any  aircraft  flying  at  a  supersonic  speed 
V  emits  a  concentrated  “boom’  -like  noise  along  rays  (Figure  .'3)  in  the  Mach  direction  (11). 
1  sketch  the  theory  of  stipersonic  booms  with  the  atmosphere  approximated  as  isothermal 
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(so  that  the  undisturbed  sound  speed  takes  a  constant  value  c  even  though  the  undisturbed 
density  p  varies  with  altitude):  a  case  permitting  quite  a  simple  extension  of  the  nonlinear 
analysis  of  waveform  shearing  and  shock  formation  (WF  and  Whitham,  1956).  Then  the  rays 
continue  as  straight  lines  at  the  Mach  angle  for  reasons  summarized  in  the  explanatory  note 
below  expression  (11).  (Actually,  the  slight  refraction  of  rays  by  temperature  stratification 
in  the  atmosphere,  when  taken  into  account  in  a  generalized  version  of  the  theory,  produces 
only  somewhat  minor  modifications  of  the  results.) 

As  such  straight  rays  stretch  out  from  a  straight  flight  path  along  cones  with  semi-angle 
(11),  any  narrow  tube  of  rays  has  its  cross-sectional  area  A  increasing  in  proportion  to 
distance  r  along  the  tube  (WF  and  Lighthill,  1956).  On  linear  theory  (WF),  acoustic  energy 
flux  u^pcA  is  propagated  unchanged  along  such  a  ray  tube  (so  that  u{pryl^  is  unchanged) 
where  u  is  air  velocity  along  it.  On  nonlinear  theory,  u{pryl^  is  propagated  unchanged  but 
at  a  signal  speed  altered  to 

(14) 


by  self-convection  and  excess- wavespeed  effects. 

This  property  can  be  described  (WF)  by  an  equation 


r  1  7  -I-  1  u  I  ^  ^1  J ,,  ... 


(WF,  pp.  187-)  where  the  quantity  in  braces  is  the  altered  value  of  the  reciprocal  of  the 
signal  speed  (14).  Now  a  simple  transformation  of  variables 


X\  =  r  —  ct^  t\  =  f  (pr)  u\  =  ^  ^  ^  — (pr)^^^  (16) 

Jo  2  c 

is  f  )und  to  convert  Equation  (15)  into  the  familiar  form 

dux  ,dux 

^  +  =  0  (17) 

atx  dxx 

which  describes  the  waveform  shearing  at  a  uniform  rate  that  is  associated  with  shock  for¬ 
mation  and  propagation  in  nonlinear  plane-wave  acoustics  (Figure  5). 

From  amongst  this  equation’s  physically  relevant  solutions  -  namely,  those  with  area- 
conserving  discontinuities  (representing  shocks)  -  the  famous  A^-wave  solution  is  the  one 
produced  by  an  initial  signal  (such  as  an  aircraft’s  pcissage  through  the  air)  that  is  first 
compressive  and  then  expansive.  The  rules  (WF)  governing  N-wave  solutions  of  Equation 
(17)  are  that  the  discontinuity  Aui  at  each  shock  falls  off  like  ix  '  while  the  space  (change 
Ax]  in  xi)  between  shocks  increases  like  .  These  rules  for  the  transformed  variables  (16) 
have  the  following  consequences  for  the  true  physical  variables:  at  a  large  distance  r  from 


the  flight  path  the  velocity  change  Au  at  each  shock  and  the  time  interval  At  between  the 
two  shocks  vary  as 


Au 


rr  *]  —1/2  r  "I 

{pr)  J  {pr)~^^^dr  and  At  k,  J  [pr)~^^^dr 


(18) 


On  horizontal  rays  (at  the  level  where  the  aircraft  is  flying),  p  is  independent  of  r  and  the 
Equations  (18)  take  the  greatly  simplified  form 


Au  «  r  and  At  « 


(19) 


appropriate  to  conical  A^-waves  in  a  homogeneous  atmosphere  (which  supersonically  con- 
vected  eddies  may  also  generate).  Actually,  the  rules  (19)  apply  also  to  the  propagation 
of  cylindrical  blast  waves  generated  by  an  exploding  wire;  since,  here  also,  ray  tube  areas 
increase  in  proportion  to  r. 

On  downward  pointing  rays  in  an  isothermal  atmosphere  p  increases  exponentially  in 
such  a  way  that  the  time  interval  At  between  shocks  approaches  the  constant  value  obtained 
in  (18)  by  making  the  integral’s  upper  limit  infinite  (WF  and  Lighthill,  1956).  On  the 
other  hand  the  shock  strength  (proportional  to  the  velocity  change  Au)  includes  the  factor 
(pr)~^/^  where  the  large  increase  in  p  from  the  flight  path  to  the  ground  (as  well  as  in 
r)  enormously  attenuates  the  supersonic  boom.  Below  Concorde  cruising  at  Mach  2,  for 
example,  an  observer  on  the  ground  hears  two  clear  shocks  with  an  interval  of  around  0.5s 
between  them,  and  yet  with  strengths  Ap/p  only  about  0.001. 


5.  Propagation  of  Sound  through  Steady  Mean  Flows 
5.1.  Adaptations  of  ray  acoustics 

Useful  information  on  sound  propagation  (including  sound  of  aeroacoustic  origin)  through 
steady  mean  flows  (Blokhintsev,  1956  and  Lighthill,  1972)  can  be  obtained  by  adaptations 
of  the  ray-acoustics  approximation.  I  sketch  these  here  before,  first,  applying  them  (in  §5.3 
below)  to  propagation  through  sheared  stratified  winds  and,  secondly,  giving  indications 
of  how  effects  of  such  parallel  mean  flows  are  modified  at  wavelengths  too  large  for  the 
applicability  of  ray  acoustics. 

Sound  propagation  through  a  steady  airflow  represents  an  autonomous  mechanical  sys¬ 
tem;  one  governed  by  laws  that  do  not  change  with  time.  Then  small  disturbances  can  be 
Fourier-analyzed  in  the  knowledge  that  propagation  of  signals  with  different  frequencies  u; 
must  proceed  without  exchange  of  energy  between  them. 
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General  Theory  of  Such  Systems  (WF,  pp.  317-) 


Such  disturbances  of  frequency  uj  involve  pressure  changes  in  the  form  P  cos  a  where  P 
varies  with  position  and  the  phase  a  is  a  function  of  position  and  time  satisfying 


da 

dt 


da 

=  u  :  the  frequency;  and  —  —  =  ki 

dxi 


the  wavenumber, 


(20) 


a  vector  with  its  direction  normal  to  crests  and  its  magnitude  27r  divided  by  a  local  wave¬ 
length. 

In  ray  theory  for  any  wave  system  (WF  and  Lighthill,  1972),  we  assume  that  the  wave¬ 
length  is  small  enough  (compared  with  distances  over  which  the  medium  -  and  its  motion, 
if  any  -  change  significantly)  for  a  well-defined  relationship 


u)  =  Q{ki,Xi) 


(21) 


to  link  frequency  with  wavenumber  at  each  position.  Equations  (20)  and  (21)  require  that 


dkj  d^a 


du  do, 


dt  dxjdt  dxj  dki 


d'^a 


+ 


d^  dil  dkj  dil 
"  + 


dxidxj )  dxj  dki  dxi  dxj  ’ 


yielding  the  basic  law  (in  Hamiltonian  form)  for  any  wave  system: 


dn 


.  ,  .  dxi  do.  .  dkj 

on  rays  satisfying  -j-  =  wavenumbers  vary  as  -7^  =  »  , 

dt  dkj  dt  dxj 


(22) 


(23) 


equations  easy  to  solve  numerically  for  given  initial  position  and  wavenumber.  However,  the 
variations  (23)  of  wavenumber  (“refraction”)  produce  no  change  of  frequency  along  rays: 


du>  _  do.  dki  dfl  dxi  dQ  f  \  dfl  dfl 

dt  dki  dt  dxi  dt  dki  \  dxi )  dxi  dki 


(24) 


so  that  rays  are  paths  of  propagation  of  the  excess  energy,  at  each  frequency,  associated  with 
the  waves’  presence. 

For  sound  waves  we  write  k  as  the  magnitude  of  the  wavenumber  vector,  expecting  that 
at  any  point  the  value  of  the  relative  frequency  in  a  frame  of  reference  moving  at  the  local 
steady  flow  velocity  u/,-  will  be  cjk  (the  local  sound  speed  times  k);  this  implies  (WF  and 
Lighthill,  1972)  that 


da 


da 


Wr  =  ^  +  Ufi—  =  U} 

dt  dxi 


Ufiki,  giving  to  =  ufr  +  ujiki  =  cjk  +  ujiki 


(25) 


as  the  acoustic  form  of  the  relationship  (21). 

[Note:  this  rule  (25)  for  relative  frequency  agrees  with  the  Doppler  rule  (6),  since  the 
velocity  of  a  source  of  frequency  u  relative  to  stationary  fluid  into  which  it  radiates  is  minus 
the  velocity  of  the  fluid  relative  to  a  frame  in  which  the  acoustic  frequency  is  u.] 
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(26) 


Use  of  this  form  (25)  of  the  relationship  (21)  in  the  ba^ic  law  (23)  tells  us  that 


dkj 

dt 


dxi 


•*u  dxi  ki 

onrayswth— =  c,-  +  u„i 


where  the  last  terms  in  these  equations  represent  adaptations  of  ray  acoustics  associated 
with  the  mean  flow.  For  example,  the  velocity  of  propagation  along  rays  is  the  vector  sum 
of  the  mean  flow  velocity  u/,  with  a  wave  velocity  of  magnitude  cj  and  direction  normal  to 
crests. 


5.2.  Energy  exchange  between  sound  waves  and  mean  flow 

The  excess  energy  (say,  E  per  unit  volume)  associated  with  the  presence  of  sound  waves 
is  propagated  along  such  rays;  in  particular,  if  attenuation  of  sound  energy  is  negligible,  then 

flux  of  excess  energy  along  a  ray  tube  =  constant.  (27) 


Note:  this  excess  energy  density  E  is  by  no  means  identical  with  the  sound  waves’  energy 
density 

Er  =<  ^-pfUsiUsi  >  +  <  ]jc)pfpl  >=  c)pf  <pI>  (28) 

(where  the  subscript  s  identifies  changes  due  to  the  sound  waves  and  the  equality  of  the 
kinetic  and  potential  energies  makes  Er  simply  twice  the  latter)  in  a  frame  of  reference 
moving  at  the  local  flow  velocity  (compare  the  definition  (25)  of  w^).  The  kinetic-energy 
part  of  the  excess  energy  density  E  is 


(29) 


which  includes  an  extra  term 


C  f 

PgUjiUji  >  =  <  pgUji  Ps~J~  ^ 

Pi  k 


(30) 


and  E  is  the  sum  of  expressions  (28)  and  (30),  giving 


U)  .  .  ,  Ur  Er  E 

E  =  Er  — ;  or,  equivalently,  either  Er  =  E —  or  —  =  — . 

Ur  U  Ur  U 


(31) 


The  quantity  E/o;,  called  action  density  in  Hamiltonian  mechanics,  is  identical  in  both 
frames  of  reference,  and  Equations  (24)  and  (27)  tell  us  that  its  flux  along  a  ray  is  constant 
(WF  and  Lighthill,  1972). 

But  Equation  (31)  shows  too  that  energy  is  exchanged  between  (i)  the  acoustic  motions 
relative  to  the  mean  flow  and  (ii)  the  mean  flow  itself.  For  example,  where  sound  waves  of 
frequency  u  enter  a  region  of  opposing  flow  (or  leave  a  region  where  the  mean  flow  is  along 
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their  direction  of  propagation)  the  ratio  uj-rlu)  increases  and  so  therefore  does  Erl E.  the 
sound  waves  gain  energy  at  the  expense  of  the  mean  flow. 

The  rate  of  exchange  of  energy  takes  the  value  (1)  written  down  in  §1.  This  is  readily 
seen  from  the  laws  governing  motion  in  an  accelerating  frame  of  reference,  which  feels 

an  inertial  force  =  -(mass)  x  (acceleration  of  frame).  (32) 

If  we  use  at  each  point  of  space  a  local  frame  of  reference  moving  with  velocity  «/,  then 
fluid  in  that  frame  has  velocity  u^i  but  is  subject  to  an  additional  force  (32);  where,  per  unit 
volume,  mass  is  pj  and  the  frame’s  acceleration  is 


dufi  ,  OUfi  .  .  ,  ^  . 

Usj-TT^  giving  force  —  PjUsj— —  doing  work  —  pf  <  UsiUsj  >  -x—  l-joj 

OXj  OXj  OXj 

per  unit  time  on  the  local  relative  motions.  This  rate  of  energy  exchange  (33)  proves  to  be 
consistent  with  the  fact  that  it  is  the  flux,  not  of  Er  but  of  action  Er/u}T,  that  is  conserved 
along  ray  tubes. 

Energy  can  be  extracted  from  a  mean  flow,  then,  not  only  by  turbulence  but  also  by 
sound  waves;  and,  in  both  cases,  the  rate  of  extraction  takes  the  same  form  (33)  in  terms  of 
perturbation  velocities  It  represents  the  effect  (§1)  of  the 

mean  momentum  flux  p/  <  Usi^sj  >  (34) 

or  Reynolds  stress  (Reynolds,  1895)  with  which  either  the  sound  waves  or  the  turbulent 
motions  act  upon  the  mean  flow.  For  sound  waves,  by  Equation  (28)  for  Er  and  by  the 

substitution  ^  ^ 

u,i  =  —p,-r,  mean  momentum  flux  =  Er^k^',  (35) 


itj,  =  —pj^i  mean  momentum  flux  =  Er-^k^', 
Pf  k  Kj 


so  that  the  Reynolds  stress  is  a  uniaxial  stress  in  the  direction  of  the  wavenumber  vector 
having  magnitude  Er- 

[Note:  strictly  speaking,  the  complete 


radiation  stress  Er 


]  for  sound  waves 

2  ‘V 


includes  not  only  the  momentum  flux  (35)  but  also  the  waves’  mean  pressure  excess 

acting  equally  in  all  directions  (Bretherton  and  Garrett,  1968);  however  (§1)  this  isotropic 
component  produces  no  energy  exchange  with  solenoidal  mean  flows.] 
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5.3.  Propagation  through  sheared  stratified  winds 

The  extremely  general  ray-acoustics  treatment  outlined  above  for  sound  propagation 
through  fluids  in  motion  has  far-reaching  applications  (in  environmental  and,  also,  in  engi¬ 
neering  acoustics)  which,  however,  are  illustrated  below  only  by  cases  of  propagation  through 
parallel  flows,  with  stratification  of  velocity  as  well  as  of  temperature  (WF  and  Lighthill, 
1972).  The  xi-direction  is  taken  as  that  of  the  mean  flow  velocity  V{xz)  which,  together 
with  the  sound  speed  c(x3),  depends  only  on  the  coordinate  X3.  Thus,  V  replaces  uj\  in  the 
general  theory  while  c  replaces  c/  (and,  for  atmospheric  propagation,  X3  is  altitude).  [Note: 
the  analysis  sketched  here  is  readily  extended  to  cases  of  winds  veering  with  altitude,  where 
Uf2  as  well  as  Uj]  is  nonzero.] 

Either  the  basic  law  (23)  or  its  ray-acoustics  form  (26)  provides,  in  general,  “refraction” 
information  in  the  form  of  three  equations  for  change  of  wavenumber;  while  the  single,  far 
simpler.  Equation  (24)  is  a  consequence  of,  but  is  by  no  means  equivalent  to,  those  three. 
By  contrast,  in  the  particular  case  when  and  cy  are  independent  of  x\  and  X2,  Equations 
(24)  and,  additionally,  (26)  in  the  cases  j  =  \  and  2  give  three  simple  results, 


u  =  constant,  /ji  =  constant  and  kz  =  constant  along  rays,  (38) 


that  may  be  shown  fully  equivalent  to  the  basic  law. 

If  now  we  write  the  wavenumber  (a  vector  normal  to  crests)  as 

{ki^  k2,  kj)  —  (k cos  V’,  «sin  ip,  k cot  0),  (39) 


so  that  K  is  its  constant  horizontal  resultant,  ip  its  constant  azimuthal  angle  to  the  wind 
direction,  and  0  its  variable  angle  to  the  vertical  (Figure  6),  and  use  Equation  (25)  in  the 
form 

U  =  c{x3)k  -t-  V{X3)ki  =  c{x3)KCOSec6  -f-  V{x3)kCOSiP,  (40) 


we  obtain  the  important 


extension  sin  6  = 


_ C(X3) _ 

U)K~^  —  V{X3)  COS  Ip 


to  Snell’s  Law 


(41) 


from  the  classical  case  (1/  =  0)  when  the  denominator  is  a  constant.  This  extended  law  (41) 
tells  us  how  9  varies  with  X3  along  any  ray  -  whose  path  we  can  then  trace,  using  Equations 
(26)  in  the  form 


=  c{x3)  cos  ips\n6  +  V  {xz),  =  c(x3)  sin  ip  sin  9, 

at  at 


dxz 

dt 


c(x3)  COS  9, 


(42) 
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by  simply  integrating  dxifdxs  and  dx2/dx3  with  respect  to  X3. 

It  follows  that  a  ray  tube  covers  the  same  horizontal  area  at  each  altitude,  so  that 
conservation  of  the  flux  of  wave  action  Erj uir  along  it  implies  that  the  vertical  component 

{Er/ojr){dx3/dt)  =  ErK~^sm6cos6  (43) 

of  wave  action  flux  is  constant  along  rays;  from  which,  with  Equation  (28),  sound  amplitudes 
are  readily  derived. 

Wind  shear  is  able  to  reproduce  all  the  main  types  of  ray  bending  (WF)  associated 
with  temperature  stratification,  and  often  to  an  enhanced  extent.  Roughly,  the  downward 
curvature  of  near-horizontal  rays  in  (km)“^  comes  to 

3  {  V'{x3)  cos  xj)  +  c\x3)}  (44) 

where  the  velocity  gradients  are  in  and  the  factor  3(km)~'s  outside  the  braces  is  an 
approximate  reciprocal  of  the  sound  speed  (WF  and  Lighthill,  1972). 

Cases  when  (44)  is  negative:  curvature  is  upward;  its  magnitude  with  zero  wind  is  at  most 
0.018  (km)-'  (because  temperature  lapse  rate  in  stable  atmospheres  cannot  exceed  10“C  per 
km,  giving  c'  =  — O.OOGs”')  but  with  strong  wind  shear  can  take  much  bigger  values  for 
upwind  propagation  (xj)  —  ir).  In  either  case  Figure  7(a)  shows  how  the  lowest  ray  emitted 
by  a  source  “lifts  off”  from  the  ground,  leaving  below  it  a  zone  of  silence  (on  ray  theory  - 
actually,  a  zone  where  amplitudes  decrease  exponentially  with  distance  below  that  ray). 

Cases  when  (44)  is  positive:  curvature  is  downward,  as  found  with  zero  wind  in  temperature- 
inversion  conditions  (e.g.,  over  a  calm  cold  lake)  and  even  more  with  strong  wind  shear  for 
downwind  propagation  {xp  =  0).  Figure  7(b)  shows  how  this  leads  to  signal  enhancement 
through  multiple-path  communication. 

In  summary,  then,  the  very  familiar  augmentation  of  sound  levels  downwind,  and  diminu¬ 
tion  upwind,  of  a  source  are  effects  of  the  wind’s  shear  (increase  with  altitude). 


5.4.  Wider  aspects  of  parallel-flow  acoustics 

The  propagation  of  sound  through  parallel  flows  at  wavelengths  too  great  for  the  appli¬ 
cability  of  ray  acoustics  can  be  analyzed  by  a  second-order  ordinary  differential  equation. 
Thus,  a  typical  Fourier  component  of  the  sound  pressure  field  takes  the  form 


{u}-Vkr)\ 


Ps  =  0. 


(45) 


Equation  (45)  can  be  used  to  improve  on  ray  acoustics 
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(a)  near  caustics  (envelopes  of  rays)  where  it  allows  a  uniformly  valid  representation  of 
amplitude  in  terms  of  the  famous  Airy  function,  giving  “beats”  between  superimposed 
waves  on  one  side  of  the  caustic,  and  exponential  decay  on  the  other  (WF); 

(b)  at  larger  wavelengths  by  abandoning  ray  theory  altogether  in  favor  of  extensive  nu¬ 
merical  solutions  of  Equation  (45);  and 

(c)  to  obtain  waveguide  modes  for  sound  propagation  in  a  two-dimensional  duct  (between 
parallel  planes)  (Pridmore-Brown,  1958;  Mungur  and  Gladwell,  1969;  and  Shankar, 
1971). 

On  the  other  hand,  in  the  case  of  a  three-dimensional  duct  carrying  parallel  flow  V{x2,  Xj) 
in  the  xi-direction,  Equation  (45)  is  converted  into  a  partial  differential  equation  (the  first 
term  being  supplemented  by  another  with  djdx^  replacing  dfdx^,  while  ^2  is  deleted)  which 
is  used 

(d)  to  obtain  waveguide  modes  in  such  ducts; 

(e)  in  calculations  of  propagation  of  sound  through  the  wide  jets  -  modelled  as  parallel 
flows  -  typical  (§4.1)  of  modern  aero-engines;  and, 

(f)  with  aeroacoustic  source  terms  included,  in  certain  enterprising  attempts  at  modelling 
jet  noise  generation  and  emission  (Phillips,  1960  and  Mani,  1976). 

6.  Acoustic  Streaming 

6.1.  Streaming  as  a  result  of  acoustic  attenuation 

Sound  waves  act  on  the  air  with  a  Reynolds  stress  (34)  even  when  mean  flow  is  absent  (so 
that  subscript  /  becomes  subscript  zero).  The  y-component  of  force  acting  on  unit  volume 
of  air  (WF,  pp.  337-)  is  then 

Q 

Ej  =  -—  <  poUsiU^j  >:  (46) 

the  force  generating  acoustic  streaming  (Lighthill,  1978a). 

However,  the  force  (46)  could  not  produce  streaming  for  unattenuated  sound  waves; 
indeed,  their  linearized  equations  can  be  used  to  show  that 

1  1  d 

if  =<  -cIp'^p]  -  -poU,iUsi  >  then  Fj  -  =<  >  (47) 
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which  is  necessarily  zero  (as  the  mean  value  of  the  rate  change  of  a  bounded  quantity). 
Accordingly,  the  fluid  must  remain  at  rest,  responding  merely  by  setting  up  the  distribution 
of  mean  pressure  whose  gradient  can  balance  the  force.  [Note:  actually,  on  the  ray- 
acoustic  approximation  (28),  is  itself  zero,  but  the  above  argument  does  not  need  to  use 
this  approximation.] 

Attenuation  of  sound  waves  takes  place 

(a)  in  the  bulk  of  the  fluid  through  the  action  of  viscosity,  thermal  conductivity  and  lags 
in  attaining  thermodynamic  equilibrium  (Ch.  6);  and 

(b)  near  solid  walls  by  viscous  attenuation  in  Stokes  boundary  layers. 

All  these  effects  produce  forces  (46)  which  act  to  generate  acoustic  streaming.  It  is 
important  to  note,  furthermore,  that  even  the  forces  due  solely  to  viscous  attenuation  - 
being  opposed  just  by  the  fluid’s  own  viscous  resistance  -  generate  mean  motions  which  do 
not  disappear  cis  the  viscosity  p  tends  to  zero  (Rayleigh,  1896;  Nyborg,  1953;  Westervelt, 
1953;  and  Nyborg,  1965). 

6.2.  Jets  generated  by  attenuated  acoustic  beams 

Attenuation  of  type  (a)  produces  a  streaming  motion  ujj  satisfying 

pufidufj/dxi  =  Fj  —  dpfdxj  -f-  pV^Ufj.  (48) 

Substantial  streaming  motions  can  be  calculated  from  this  equation  only  with  the  left-hand 
side  included  (Stuart,  1966);  although  in  pre-1966  literature  it  was  misleadingly  regarded  as 
“a  fourth-order  term”  and  so  ignored  -  thus  limiting  all  the  theories  to  uninteresting  cases 
when  the  streaming  Reynolds  number  would  be  of  order  1  or  less. 

We  can  use  streaming  generated  by  acoustic  beams  to  illustrate  the  above  principles.  If 
acoustic  energy  is  attenuated  at  a  rate  /?  per  unit  length,  then  a  source  at  the  origin  which 
beams  power  P  along  the  zi-axis  transmits  a  distribution  of 

power  ,  and  therefore  energy  per  unit  length  c~^  ;  (49) 

which  is  necessarily  the  integral  of  energy  density,  and  so  also  of  the  uniaxial  Reynold  stress 
(35),  over  the  beam’s  cross-section.  It  follows  by  differentiation  that  the  force  per  unit 
volume  (46),  integrated  over  a  cross-section,  produces  (Lighthill,  1978a) 

a  force  per  unit  length  in  the  xj-direction.  (50) 

At  high  ultrasonic  frequencies  the  force  distribution  (50)  is  rather  concentrated,  the  distance 
of  its  center  of  application  from  the  origin  being  just  (which  at  1  MHz,  for  example,  is 
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24mm  in  air).  Effectively  (WF,  p.  345),  the  beam  applies  at  this  center  a  total  force  c  ^  P 
(integral  of  the  distribution  (50)). 

The  type  (WF  and  Lighthill,  1978a)  of  streaming  motion  generated  by  this  concentrated 
force  c~'P  depends  critically  on  the  value  of  a  sort  of  Reynolds  number  squared, 

which  is  about  lO^P  in  atmospheric  air  (with  P  in  watts).  Streaming  of  the  low-Reynolds- 
number  “stokeslet”  type  predicted  (for  a  concentrated  force)  by  Equation  (48)  with  the  left- 
hand  side  suppressed  is  a  good  approximation  only  for  P  <  For  a  source  of  power 

10“'* VF,  by  contrast,  the  force  c"*P  generates  quite  a  narrow  laminar  jet  with  momentum 
transport  c“*  P,  and  at  powers  exceeding  3  x  10“'‘VF  ‘this  jet  has  become  turbulent,  spreading 
conically  with  semi-angle  about  15°  and  continuing  to  transport  momentum  at  the  rate  c“’P. 
Such  turbulent  jets  gen^^rated  by  sound  are  strikingly  reciprocal  to  a  classical  aeroacoustic 
theme! 

At  lower  frequencies  an  acoustic  beam  of  substantial  power  delivers  a  turbulent  jet  with 
a  somewhat  more  variable  angle  of  spread  -  but  one  which 

at  each  point  x\  carries  momentum  transport  c“*P(l  —  e"^^‘),  (51) 

generated  by  the  total  force  (50)  acting  up  to  that  point.  This  momentum  transport  in  the 
jet  represents  the  source’s  original  rate  of  momentum  delivery  minus  the  acoustic  beam’s 
own  ren^aining  momentum  transport  (49).  In  summary,  as  acoustic  power  is  dissipated  into 
heat,  the  associated  acoustic  momentum  transport  is  converted  into  a  mean  motion  (which, 
at  higher  Reynolds  numbers,  is  turbulent:  Lighthill,  1978a). 


6.3.  Streaming  around  bodies  generated  by  boundary-layer  attenuation 


Sound  waves  of  frequency  w  well  below  high  ultrasonic  frequencies  have  their  attenuation 
concentrated,  if  solid  bodies  are  present,  in  thin  Stokes  boundary  layers  attached  to  each 
body  (WF).  Then  the  streaming  generated  near  a  particular  point  on  a  body  surface  is  rather 
simply  expressed  by  using  local  coordinates  with  that  point  as  origin,  with  the  2-axis  normal 
to  the  body  and  the  i-axis  in  the  direction  of  the  inviscid  flow  just  outside  the  boundary 
layer:  the  exterior  flow.  The  Stokes  boundary  layer  for  an  exterior  flow 


{U{x,y),V{x,y))e*'^^  has  interior  flow 
{U(x,y),V{x,  [l  -  . 


(52) 


[Note  that  my  choice  of  coordinates  makes  F(0,0)  =  0,  and  that  the  expressions  (52)  become 
identical  outside  the  layer.]  The  streaming  (WF  and  Lighthill,  1978a)  is  calculated  from  the 
equation 


P/nt  _  fExr  ^  0^ 


(53) 
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with  certain  differences  from  Equation  (48)  explained  as  follows: 

(a)  the  first  term  is  the  force  (46)  generating  streaming  within  the  boundary  layer; 

(b)  we  are  free,  however,  to  subtract  the  second,  since  (see  §6.1)  it  can  produce  no  stream¬ 
ing,  and  conveniently,  the  difference  is  zero  outside  the  layer; 

(c)  gradients  in  the  ^-direction  are  so  steep  that  the  third  term  dominates  the  viscous  force 
-  and,  indeed,  in  such  a  boundary  layer,  dominates  also  the  left-hand  side  of  Equation 
(48). 

The  solution  of  Equation  (53)  which  vanishes  a,t  z  —  0  and  has  zero  gradient  at  the  edge 
of  the  layer  is  obtained  by  two  integrations,  and  its  exterior  value  is 

“"T  =  f‘-'  r  (fr  -  (54) 

where  integration  extends  in  practice,  not  to  “infinity,”  but  to  the  edge  of  the  layer  within 
which  the  integrand  is  nonzero.  Expression  (54)  for  the  exterior  streaming  is  yet  again  (see 
§6.1)  independent  of  the  viscosity  since  Equation  (52)  makes  zdz  of  order  n/ pu\  and  it  is 
easily  evaluated. 

At  I  =  ?/  =  0  (in  the  coordinates  specified  earlier)  the  exterior  streaming  (54)  has 

^  ,,3dUldx  +  2dV/dy  ,  ^  ,^dV/dx 

X  —  component  —  U - - - and  t/— component  —  U - ,  (55) 

4u;  Au) 

with  zero  2-component.  This  is  a  generalized  form  of  the  century-old  Rayleigh  law  of  stream¬ 
ing  (which  covers  cases  when  V  is  identically  zero). 

For  the  complete  streaming  pattern,  expressions  (55)  are,  effectively,  boundary  values 
for  its  tangential  component  at  the  body  surface  (because  the  Stokes  boundary  layer  is  so 
thin).  Therefore,  any  simple  solver  for  the  steady-flow  Navier-Stokes  equations  with  specified 
tangential  velocities  on  the  boundary  allows  the  pattern  to  be  determined.  Important  note: 
here,  the  inertial  terms  in  the  Navier-Stokes  equations  must  not  be  neglected,  unless  the 
Reynolds  number  /?,  ba.sed  on  the  streaming  velocity  (55)  be  of  order  1  or  less;  when, 
however,  the  corresponding  streaming  motions  would  (as  in  §6.2)  be  uninterestingly  small. 

In  the  other  extreme  case  when  /?,  is  rather  large  (at  least  10^)  the  streaming  motion 
remains  quite  close  to  the  body  (Stuart,  1966)  within  a  steady  boundary  layer  whose  di¬ 
mension  (relative  to  that  of  the  body)  is  of  order  This  layer  is  by  no  means  as  thin 

as  the  Stokes  boundary  layer,  but  it  does  confine  very  considerably  the  acoustic  streaming 
motion.  Equations  (55)  direct  this  motion  towards  one  of  the  exterior  flow’s  stagnation 
points,  whence  the  steady-boundary-layer  flow  emerges  as  a  jet  (Figure  8)  -  yet  another  jet 
generated  by  sound  (Riley,  1987  and  Amin  and  Riley,  1990). 
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Figure  1 


(a)  When  a  source  of  sound  of  period  T  ap-  (b)  Case  of  observer  at  angle  6  to  source’s 

proaches  an  observer  at  velocity  w,  succes-  direction  of  motion  at  speed  V ;  then  veloc- 

sive  wave  crests  are  emitted  towards  him  ity  of  approach  is  u;  =  V  cos  d. 

with  separation  \  =  cT  —  wT. 


Figure  2.  Sounds  arriving  simultaneously  from  a  source’s  far  side  F  and  near  side  N  took 
times  t  and  t  —  t  (say)  to  reach  the  observer,  while  the  distance  FN  increased  from  (  to 
i  +  WT  during  the  time  r  between  emissions  from  F  and  N. 
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Figure  4.  A  uniformly  valid  Doppler-effect  approximation. 

Diagram  (a)  Space-time  diagram  for  unconvected  “eddies”  of  correlation  length  i  and 
duration 

Diagram  (b)  Ccise  of  “eddies”  convected  towards  observer  at  velocity  w\  being  Diagram 
(a)  sheared  by  a  distance  w  per  unit  time.  Here,  lines  sloping  by  a  distance  c  per  unit  time 
represent  emissions  received  simultaneously  by  observer. 

Case  (i):  wjc  small.  Case  (ii):  wjc  =  1.  Case  (iii):  intermediate  value  of  u;/c. 

Diagram  (c) - Average  modification  factor  (bl)- 

. Acoustic  efficiency,  obtained  by  applying  this  factor  to  a  low-Mach-number 

“quadrupole”  efficiency  of  (say)  10~^(l//c)^. 
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Figure  5.  Equation  (17)  tells  us  that  any  value  of  ui  propagates  unchanged  along  a  char¬ 
acteristic  dxi/dti  =  ui.  This  implies  (WF,  p.  151)  waveform  shearing  at  a  uniform  rate 
(Diagram  (a)).  According  to  nonlinear  acoustics  (WF,  pp.  170-)  area-conserving  discontinu¬ 
ities,  shown  as  broken  lines  in  Diagram  (a),  have  to  be  incorporated  wherever  necessary  to 
keep  the  solution  one-valued;  leading,  in  the  case  illustrated,  to  the  N-w&ve  form  (Diagram 

(b)). 


I 


Figure  6.  The  angle  6  between  a  ray  and  the  vertical  varies  in  accordance  with  the  extension 
(41)  to  Snell’s  law. 


(a)  (b) 


Figure  7.  Effects  of  ray  curvature  (44)  on  propagation  from  a  source  on  horizontal  ground. 

Diagram  (a)  Rays  of  given  upward  curvature  (due  to  temperature  lapse  or  upwind  prop¬ 
agation)  can  leave  a  Zone  of  Silence  (ZS)  below  the  ray  emitted  horizontally. 

Diagram  (b)  Rays  of  given  downward  curvature  (due  to  temperature  inversion  or  down¬ 
wind  propagation)  can  enhance  received  signals  through  multiple-path  communication. 
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Incident 
sound  waves 


Figure  8.  Plane  sound  waves  incident  on  a  sphere  may  generate  a  steady  streaming  motion 
(Amin  and  Riley,  1990),  concentrated  in  a  relatively  thick  boundary  layer,  and  directed 
towards  one  of  the  exterior  flow’s  stagnation  points  -  whence  it  emerges  as  a  jet. 
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